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ABSTRACT 

Weak lensing convergence can be used directly to map and probe the dark mass distri- 
bution in the universe. Building on earlier studies, we recall how the statistics of the 
convergence field are related to the statistics of the underlying mass distribution, in 
particular to the many-body density correlations. We describe two model- independent 
approximations which provide two simple methods to compute the probability distri- 
bution function, pdf, of the convergence. We apply one of these to the case where the 
density field can be described by a log- normal pdf. Next, we discuss two hierarchical 
models for the high-order correlations which allow one to perform exact calculations 
and evaluate the previous approximations in such specific cases. Finally, we apply 
these methods to a very simple model for the evolution of the density field from lin- 
ear to highly non-linear scales. Comparisons with the results obtained from numerical 
simulations, obtained from a number of different realizations, show excellent agree- 
ment with our theoretical predictions. We have probed various angular scales in the 
numerical work and considered sources at 14 different redshifts in each of two dif- 
ferent cosmological scenarios, an open cosmology and a flat cosmology with non-zero 
cosmological constant. Our simulation technique employs computations of the full 3-d 
shear matrices along the line of sight from the source redshift to the observer and is 
complementary to more popular ray-tracing algorithms. Our results therefore provide 
a valuable cross-check for such complementary simulation techniques, as well as for 
our simple analytical model, from the linear to the highly non-linear regime. 

Key words: Cosmology: theory - gravitational lensing - large-scale structure of 
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1 INTRODUCTION 

Weak gravitational lensing is responsible for the shearing 
and magnification in the images of high-redshift sources due 
to the presence of intervening mass. Since the lensing effects 
arise from deflections of the light rays due to fluctuations of 
the gravitational potential, they can be directly related to 
the underlying density field of the large-scale structure and 
not necessarily to the presence of luminous structure in the 
form of galaxies and clusters. Consequently, statistical anal- 
ysis of observed weak lensing data has proved to be a pow- 
erful tool to probe the underlying density, which is assumed 
to be dominated by the dark matter content. Since the un- 
derlying density field and its evolution is strongly dependent 
on the cosmological parameters, observational surveys have 
been particularly fruitful in providing estimates of the pa- 



rameters, which can be pre-set in analytical calculations or 
in numerical simulations. 

Analytical computations for weak lensing statistics can 
be readily obtained for large smoothing angles, > 10', where 
perturbative calculations are applicable (e.g., Villumsen, 

1996, Stebbins, 1996, Bernardeau et al., 1997, Jain & Seljak, 

1997, Kaiser, 1998, Van Waerbeke, Bernardeau & Mellier, 
1999, and Schneider et al., 1998). However, on small angular 
scales, < 5', especially relevant to observational surveys with 
small sky coverage, perturbative calculations are no longer 
valid and models to represent the gravitational clustering in 
the non-linear regime have had to be devised. 

Thus, a number of useful techniques have been devel- 
oped to transform from the linear matter power spectrum 
to a fully non-linear spectrum. Originally, Hamilton et al. 
(1991) considered the evolution of the matter correlation 
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function and this work was extended by Peacock & Dodds 
(1996) to describe the non- linear evolution of the matter 
power spectrum. Their methods were based on the conser- 
vation of mass and a rescaling of physical lengths in the 
different regimes, following the "stable-clustering" Ansatz 
which assumes that small scales decouple and are statisti- 
cally frozen in proper coordinates (Peebles, 1980). 

In a different approach, Peacock & Smith (2000) and 
Seljak (2000) developed a model based on the random dis- 
tribution of dark matter haloes, modulated by the large- 
scale matter distribution. This "Halo Model" for non-linear 
evolution is able to reproduce the matter power spectrum 
of iV-body simulations over a wide range of scales and has 
the advantage of relating the linear and non-linear power at 
the same scale through fitting formula?. Smith et al. (2002) 
have presented a new set of fitting functions based on the 
Halo Model and calibrated to a set of iV-body simulations. 
Barber & Taylor (2002) have shown excellent agreement be- 
tween the power spectrum in the lensing convergence ob- 
tained from numerical simulations and the predictions from 
the Halo model fitting functions of Smith et al. (2002). 

An advantage of the halo model is that in principle 
it can predict the higher-order correlation functions in the 
highly non-linear regime. Indeed, on small scales the latter 
are set by the density profile of the halos as the p— point cor- 
relation is dominated by the contribution associated with all 
p points being within the same halo. However, the neglect 
of substructures may lead to larger inaccuracies for higher- 
order statistics. Moreover, at intermediate scales one also 
probes the correlations among different halos which intro- 
duces new unknowns and makes explicit calculations cum- 
bersome. Finally, such a model for the density field is not 
well-suited to describe low-density and underdense regions 
(e.g., voids, filaments) which are outside virialized objects. 

These problems have motivated the recourse to an al- 
ternative approach which is directly based on the many- 
body correlations. The most common model of this kind 
is to express the p— point correlations as a sum of p — 1 
products over the two-point correlations linking all p points. 
This yields the class of "Hierarchical models," which are 
specified by the weights given to any such topology associ- 
ated with thep— 1 products (e.g., Fry 1984, Schaeffer 1984, 
Bernardeau & Schaeffer 1992, and Szapudi & Szalay 1993, 
1997). Once these weights have been assigned it is possible 
to resum all many-body correlations and to compute the 
probability distribution function (pdf) of the density field, 
or of any quantity which is linearly dependent on the matter 
density. This is most easily done for "minimal tree-models" , 
where the weight associated to a given tree-topology is set by 
its vertices (e.g., Bernardeau & Schaeffer 1992), or for "stel- 
lar models", which only contain stellar diagrams (Valageas, 
Barber & Munshi, 2003). Such an approach is very well- 
suited to the study of weak-lensing which only involves the 
matter density field and does not make the distinction be- 
tween astrophysical objects like clusters, Lyman-a clouds or 
voids. Hence the many-body correlations of the density field 
are precisely the quantities of interest which directly appear 
through the statistics of weak-lensing effects, rather than 
the possible decomposition of the density field over different 
classes of objects. 

Using such an approach, coupled to the Hamilton et al. 
(1991) prescription for the two-point correlation, Valageas 



(2000a, b) and Munshi & Jain (2000 and 2001) were able to 
compute the pdf of the weak-lensing convergence whilst the 
associated bias was considered by Munshi (2000) and the cu- 
mulant correlators associated with such distributions were 
evaluated by Munshi & Jain (2000). Munshi & Wang (2003) 
further extended these studies to cosmological scenarios in- 
cluding dark energy. These methods can also handle more 
intricate quantities like the aperture-mass or the shear which 
involve compensated filters and require a detailed model for 
the many-body correlations. Thus, using a minimal tree- 
model for the non-linear regime, Bernardeau and Valageas 
(2000) were able to predict the pdf of the aperture-mass 
and to obtain a good agreement with numerical simulations 
(they also showed that similar techniques could be applied 
to the quasi-linear regime where the calculations can actu- 
ally be made rigorous) . On the other hand, adopting a stellar 
model for the many-body correlations, Valageas, Barber & 
Munshi (2003) obtained excellent agreement for the shear 
pdf when compared with the results of iV-body simulations. 

A clear alternative to the analytical approaches for the 
understanding of weak lensing has been provided by the de- 
velopment of numerical techniques to simulate different cos- 
mological scenarios. A recent initiative, designed specifically 
for weak lensing studies in numerical simulations, was devel- 
oped by Couchman, Barber & Thomas (1999) whose tech- 
nique allows the computation of the full 3-dimensional shear 
matrices at locations along the lines of sight entirely within 
each simulation volume. These intermediate matrices are 
then combined using the full implementation of the multi- 
ple lens-plane theory (described fully by Schneider, Ehlers & 
Falco, 1992) to obtain the final Jacobian matrix for each line 
of sight, appropriate for sources at the selected redshift. Ex- 
cellent agreement between the results obtained using Couch- 
man et al.'s (1999) method and analytical predictions have 
been reported. Barber (2002) has compared the redshift and 
scale dependence of the shear variance with the analytical 
program of Van Waerbeke et al. (2001), based on fitting for- 
mulae. Barber & Taylor (2002) computed the angular power 
spectra for the convergence and magnification and obtained 
higher-order moments for the convergence; their results were 
in close agreement with the predictions from a halo model, 
using the functional fitting devised by Smith et al. (2002). 
Most recently, Valageas, Barber & Munshi (2003) were able 
to confirm the predictions of their hierarchical Ansatz, using 
a stellar model, for the full pdf of the shear for sources at 
redshift 1 and on angular scales representative of the non- 
linear regime. 

The main goals of the current paper are to test against 
numerical simulations a very simple model for the evolution 
of the density field which allows one to compute the statistics 
of the weak-lensing convergence from linear to highly non- 
linear scales and to compare the various approximations one 
can use within this framework. To achieve this we present re- 
sults for a selection of source redshifts and for angular scales 
from l'.O to 8'.0, taking us from the quasi-linear regime into 
the highly non-linear regime. We concentrate here on the 
full pdf and low-order moments for the lensing convergence. 
We present the following analytical approximations. 

(i) We approximate the mean of the many-body correla- 
tions over the cone, represented by the angular window, by 
the mean evaluated over spherical shells. This approxima- 
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tion is as described by Valageas (2000b) and it is completely 
model-independent. Indeed, it can be applied to any descrip- 
tion used for the density field. 

(ii) Adding a further approximation to the previous method, 
we can approximate the generating function, <p/i s (y), for the 
smoothed normalised convergence, k a , (from which the pdf is 
obtained after performing a Laplace transform) by the gen- 
crating function, tp{y), associated with the density contrast, 
5r, at the typical scale and redshift probed by observations. 
As shown in Valageas (2000b), this procedure is sufficient to 
recover the properties of the convergence with a reasonable 
accuracy and it is again model-independent. 

(iii) We apply the previous approximation to the specific 
case where the pdf of the density contrast can be described 
by a log-normal law. 

(iv) We briefly recall how to compute exactly the pdf of 
the convergence in the case where the connected correlations 
of the density field are given by a minimal-tree model, fol- 
lowing Valageas (2000b) or Bernardeau & Valageas (2000). 

(v) We show how to compute the statistics of the conver- 
gence from the stellar model for the connected correlations 
of the density field introduced in Valageas, Barber & Munshi 
(2003) to derive the pdf of the weak-lensing shear. 

Eventually we apply these methods to a very simple model 
for the evolution of the density field from the linear to the 
highly non-linear regime. 

Finally, the analytical predictions are compared with 
the results of numerical simulations, in which Couchman et 
al.'s (1999) fully 3-dimensional shear method is used. The 
numerical data have been obtained in two different cosmolo- 
gies which we will describe as LCDM (a flat cosmology with 
a cosmological constant) and OCDM (an open cosmology 
with zero cosmological constant). The data have been com- 
puted for sources at 14 different redshifts and on a wide 
range of angular scales. 

This paper is organized as follows. Section 2 defines our 
notations and outlines the basic equations which express the 
weak-lensing convergence in terms of the density field along 
the line of sight. Section 3 provides a short description of our 
analytical model for the evolution of the density field from 
linear to highly non-linear scales. We also recall the general 
relationship between the pdf of the density contrast and its 
cumulant generating function and we briefly discuss two hi- 
erarchical models which can be used to describe the many- 
body density correlations. Section 4 provides the main ana- 
lytical results and highlights various approximations which 
can be used to derive the statistics of convergence maps. 
Section 5 describes our simulations and section 6 contains 
a detailed comparison of simulation results with analytical 
predictions. Finally section 7 is reserved for discussion of our 
results within the context of observational programme and 
future plans. 



2 THE CONVERGENCE FIELD 

As is well-known, the image of a distant source at a redshift 
z a received by an observer at redshift z = is distorted by 
the deflection of light due to density fluctuations along the 
line of sight. This shearing and magnification are described 
in terms of angular position vector, by the shear tensor 
dM/dti, whose trace, n, (also called the convergence) yields 



the magnification of the source (in the weak-lensing regime) . 
The convergence is given by (Bernardeau et al., 1997, and 
Kaiser, 1998): 



3o» r 

2 Jo 



(1) 



In this equation, f2 m is the matter density parameter, x an d 
\s are the radial distance along the line of sight and radial 
distance to the source, respectively, <5(x) = (p(x) — p)/p is 
the density contrast at position x and 



Hi V(x)V( X s - x) 



(1+*) 



(2) 



ensures the integration takes account of the relevant angular 
diameter distances T>(x) to radial distance x, 2?(Xs) to the 
radial position of the source and T>(xs — x) from x to Xs! 



dz 



^/n A + (1 - On - A )(1 + «) 2 + fim(l + Z) 3 ' 



(3) 



where c is the velocity of light, Ho is the present-day value of 
the Hubble parameter and Oa is the vacuum energy density 
parameter; finally, 



v(x) 



-f- ShlK (| 1 - On - «A | 1/2 H X/C) 



(4) 



in which sin^ means the hyperbolic sine, sinh, if (1 — £3 m — 
Qa) > 0, or sine if (1 - fi m - Q A ) < 0; if (1 - fi m - fi A ) = 0, 
then T>(x) = X- Eq.(l) assumes that the components of the 
shear tensor are small so that we can use the Born approx- 
imation (i.e. the integral over redshift is taken along the 
unperturbed line of sight) but the density fluctuations S 
can be large (Kaiser, 1992). Next, we can see from eq.(l) 
that there is a ininircium value, K m i n (z s ), for the conver- 
gence of a source located at redshift z s , which corresponds 
to an "empty" beam between the source and the observer 
(5 = — 1 everywhere along the line of sight): 



3fi n 



dx w(x,Xs)- 



(5) 



Following Valageas (2000a, b) and Munshi & Jain (2000) it 
is convenient to define the "normalised" convergence, ft, by: 



dx w S, with w 



w(x,Xs) 



J Xs dx w(x,Xs) 



(6) 



which obeys ft > — 1. Here we introduced the "normalised 
selection function," w(x,Xs)- As shown in Valageas (2000a, 
b), one interest of working with normalised quantities like 
k is that most of the cosmological dependence (on f2 m ,OA 
and Zs) and the projection effects are encapsulated within 
Kmin, while the statistics of ft (e.g., its pdf) mainly probe the 
deviations from Gaussianity of the density field which arise 
from the non-linear dynamics of gravitational clustering as 
well as the amplitude of the density fluctuations (erg). If one 
smoothes the observations with a top-hat window in real 
space of small angular radius, 6 S , one rather considers the 
filtered normalised convergence k s (where the subscript "s" 
refers to "smoothed" ) : 



Jo ^ 2 Jo 



dx w{x,Xs) S (x,T>3) ■ 



(7) 
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Here $ is a vector in the plane perpendicular to the line of 
sight (we restrict ourselves to small angular windows) over 
which we integrate within the disk |$| < 8 S ; we note this by 
the short notation J Q . Thus x is the radial coordinate, while 
T>{) is the two-dimensional vector of transverse coordinates. 
Eq.(7) clearly shows that the smoothed convergence k s is 
actually an average of the density contrast over the cone of 
angular radius 6 S . 

For some purposes it is convenient to work in Fourier 
space. Thus, we define the Fourier transform of the density 
contrast by: 



dk e 1 



<5(k) 



(8) 



where x and k are comoving coordinates. Then, eq.(7) also 
reads: 



(9) 



«. = J d X w(x,X.) J dk e ife " x W(k ± V6 s ) <5(k), 

where k« is the component of k parallel to the line of sight 
and kx is the two-dimensional vector formed by the compo- 
nents of k perpendicular to the line of sight. Here we intro- 
duced the Fourier form W(k±T>8 s ) of the real-space top-hat 
filter of angular radius 9 S : 

where Ji is the Bessel function of the first kind of order 1. 
If we choose another filter (e.g., a Gaussian window rather 
than a top-hat) the expression (9) remains valid and we 
simply need to use the relevant Fourier window W(k±T>6 s ). 
In Fourier space, we also define the power-spectrum, P(k), 
of the density contrast by: 



<<y(ki)<y(ka)) = fo(ki+k 2 ) P(ki), 

where 8b is Dirac's distribution. Then, we obtain: 



(11) 



£„(*) = (J(xi)J(xi + x)> = J dk e ikx P(k), (12) 
for the two-point correlation £2(2) of the density contrast. 



3 ANALYTICAL DESCRIPTION OF THE 
DENSITY FIELD 

In order to derive the pdf of the convergence, V(k), we 
clearly need to specify the properties of the underlying den- 
sity field. As seen from eq.(l), the convergence, k, is given by 
a linear integral along the line of sight over the density field. 
Then, the direct computation of the pdf of such a sum over 
the redshift of the lenses would yield an infinite number of 
convolution products which makes it intractable. However, 
this problem can be greatly simplified by working with the 
logarithm of the Laplace transform of the pdf, <p(y). Indeed, 
since the Laplace transform changes convolutions into or- 
dinary products and the logarithm changes products into 
sums, the generating functions ifi(y) simply add when differ- 
ent layers along the line of sight are superposed. This is the 
basis of the method introduced in Valageas (2000a, b). This 
approach has already been presented in detail in various 
works, for the convergence (e.g., Valageas, 2000a, b; Munshi 



& Jain, 2000 and 2001), the aperture mass (Bernardeau & 
Valageas, 2000) and the shear (Valageas, Barber & Munshi, 
2003), using various models for the density field which are 
well-suited to such a technique. 

Therefore, we briefly review in the following sub-sections 
the various phenomenological models which have been put 
forward to describe the density field, which we aim to com- 
pare in this article through their implications for the pdf of 
the convergence, V(k). 



3.1 Cumulant generating function for the density 
contrast 

As recalled above, in order to handle the projection effects 
associated with the integration of the density fluctuations 
along the line of sight, it is convenient to work with the 
logarithm <p(y) of the Laplace transform of the pdf. There- 
fore, it is useful to define also the pdf V(Sr) of the density 
contrast at scale R through its generating function <p(y) : 



-1 



d8 R e 



(13) 



where Sr is the density contrast within spherical cells of 
radius R and volume V while £ 2 is its variance: 



5«= / 
Jv 



dx 

V 



<5(x) and Z a = {6%). 



(14) 



The pdf V(Sr) can be recovered from tp(y) through the in- 
verse Laplace transform: 



V(Sr) = 



(15) 



As is well-known, the function (p(y) defined from eq.(13) is 
also the generating function of the cumulants of the density 
contrast (see any textbook on probability theory). Thus, the 
expansion of <p{y) at y = reads: 



<p(y) = 



(-iy 



pi 



5 P r 



with S v = 



(Sr}c 

ir 



(16) 



The reason it is useful to introduce the variance £ 2 in the 
definition (13) of the cumulant generating function ip(y) is 
that it removes most of the dependence on scale and time 
of the properties of the density field. More precisely, it can 
be shown that <p(y) as defined above has a finite limit in the 
limit £ 2 — > 0, which corresponds to the quasi-linear regime. 
This exact result can be obtained from the expansion (16) 
through a perturbative method (Bernardeau, 1994) or more 
rigorously from eq.(13) through a steepest-descent method 
(Valageas, 2002). In particular, the derivation of the generat- 
ing function ip(y) in this quasi-linear limit yields the implicit 
system: 



<p(v) 



C(r)- 

-vC'(r) 



rC'(r) 



(17) 
(18) 



where the function £(r) is closely related to the spherical 
dynamics for the non-linear density contrast (Bernardeau, 
1994; Valageas, 2002). This function only depends on 
the local slope n of the power-spectrum of the density fluc- 
tuations and on the cosmological parameters il nl (z),flA(z). 
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However, the dependence on Q m , Qa, is rather small (Bernardeau, 
1992) so that over the whole quasi-linear regime (£ 2 1) 
the pdf V(5r) can be fully described through two quanti- 
ties only: the variance £2 an d the local slope n (which yields 
£(t), whence f(y) and finally V(5r) using £ 2 )- 

In the non-linear regime there are no more rigorous re- 
sults for the behaviour of the pdf P(5r). However, a rea- 
sonable approximation is provided by the "stable-clustering 

Ansatf (e.g., Peebles 1980; Balian & Schaeffer 1989) where 

— p— 1 

the cumulants obey the scaling law (8 R ) C oc £ 2 • Therefore, 
the coefficients S p introduced in eq.(16) are again constant 
with time so that the generating function <p(y) is again inde- 
pendent of £ 2 (i-6- it is unique for the whole highly non- linear 
regime £ 2 ^> 1). However, it should still depend on the local 
slope n of the power-spectrum. Note that deviations from 
the stable-clustering Ansatz simply mean that (p(y) should 
still exhibit a weak dependence on time within the highly 
non- linear regime. 



(a) 



(b) 



(c) 



3.2 A simple parameterization for the generating 
function <p(y) 

From the properties of the generating function (p(y) re- 
called above in Sect. 3.1, we choose the following sim- 
ple parameterization for ip(y). For all regimes, from quasi- 
linear scales down to highly non-linear scales, we parameter- 
ize <fi(y) through the associated function £(r) which obeys 
the implicit system (17)-(18). Next, we choose a simple 
phenomenological prescription for CM- Following previous 
works (e.g., Bernardeau & Schaeffer 1992, Bernardeau & 
Valageas 2000, Valageas, Barber & Munshi, 2003) we use: 



CM 



(19) 



where we have kept the usual notation, k (not to be confused 
with the convergence), for the free parameter which enters 
the definition of (,(t). In order to handle the variation of 
f>{y) (hence of (,{t)) from linear to highly non-linear scales, 
we let this parameter k vary so as to recover the correct 
skewness S3 of the density contrast. From eqs.(17)-(18) this 
yields: 

K = sh- ™ 

For S3 > 3 we have k > while for < S3 < 3 we have 
k < —1. In the limit S3 — * 3 the function £(t) defined in 
eq.(19) goes to £(t) — > e~ T — 1 (and \n\ —> 00). Then, we 
take for the skewness in the highly non- linear regime, S^ L , 
the prediction of HEPT (Scoccimarro & Frieman, 1999) and 
for the quasi- linear regime, S^ L , the exact result obtained 
from perturbative theory: 

2" 



sr = 3- 



5? L 



34 



(n + 3). 



(21) 



1 + 2™+!' ~ a 7 

Here n is the local slope of the linear power-spectrum at 
the typical wavenumber k s probed by the observations at 
redshift z. From eq.(9), we define this typical wavenumber 



k s (z) = 



V(z)6 s 



(22) 



Finally, we introduce the power, A 2 (fc), per logarithmic in- 
terval defined as: 





Figure 1. The various topologies one can build for the 5-point 
connected correlation function within the framework of a tree- 
model as in (27). Graph (a) is a "stellar diagram" (i.e., (p — 1) 
points arc linked to a central point), while graph (b) is a "snake 
diagram" (i.e., one successively visits all points along one curve 
with only two end-points) and graph (c) is a "hybrid diagram" . 



A 2 (k,z) = 4nk 3 P(k,z). (23) 

Then, for intermediate regimes defined by 1 < A 2 (k s ,z) < 
A v i r (z), we use the simple linear interpolation: 

S>{z) = S? L + ^ggy^ (S 3 NL - . (24) 

For A 2 (fc s , z) < 1 (i.e. the quasi-linear regime), we take S3 = 
S^ L , while for A 2 (fc s ,z) > A v i r (,z) (i.e., the highly non- 
linear regime), we take S3 = S^ L . Here A v ir is the density 
contrast at virialization given by the usual spherical collapse 
(thus A v i r ~ 178 for a critical-density universe). Indeed, the 
threshold A 2 > A v i r describes the highly non-linear regime 
where most of the matter at scale 1/fc has collapsed into 
non-linear structures. 

Note that the function £(r) which we use as a mere 
intermediate tool to parameterize the generating f{y) also 
has a more physical meaning. As recalled in Sect. 3.1, in the 
quasi-linear regime it can actually be derived from the equa- 
tions of motion and it is closely related to the spherical dy- 
namics. On the other hand, in the highly non-linear regime 
it can also be interpreted as an approximation to the vertex 
generating function £ v (t) which appears within the frame- 
work of "minimal tree-models" , as we shall discuss below in 
Sect. 3.3. Since we use a constant skewness S^ L in the highly 
non-linear regime this parameterization for tp(y) is consis- 
tent with the stable-clustering Ansatz. However, if needed it 
is straightforward to include deviations from this Ansatz by 
incorporating some additional dependence on time into the 
skewness S3. 

Finally, we obtain the non-linear evolution of the power- 
spectrum from the prescription given by Peacock & Dodds 
(1996). This completes the description of the pdf V(5r) and 
of the cumulants (<5^} c for all scales and times and for any 
cosmological parameters. 



3.3 Minimal tree model 

The generating function ip(y) introduced in Sect. 3.1 de- 
fines the properties of the density contrast 8r smoothed 
over spherical cells of radius R. In order to fully describe 
the density field we actually need to specify the detailed be- 
haviour of the many-body connected correlation functions, 
f p (xi, .., x p ), defined by (Peebles, 1980): 



£p(xi,..,x p ) = (<5(xi)..<5(x p )) c 



(25) 
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Indeed, the cumulants (5 R ) C introduced in the previous sec- 
tion only measure the mean of these correlations over spher- 
ical cells: 



P > 2 : {5 R ) C = t 



■-I 



dxi..dx p 



Cp( X ll "I X p)' 



(26) 



Therefore, the cumulants (<5^} c are insufficient if we are in- 
terested in real-space niters which are significantly different 
from a spherical top-hat. As shown in Valagcas (2000a, b), 
this is not very important for the convergence, k, but it 
is a key point for the aperture-mass, M ap , see Bernardeau 
& Valageas (2000), or the shear, 7, see Valageas, Barber 
& Munshi (2003), which involve compensated filters. This 
leads one to introduce more precise models for the density 
field which fully describe the correlations £ p (xi, .., x p ). 

One such model is the "minimal tree-model" which is 
actually a specific case of the more general "tree-models". 
The latter are defined by the hierarchical property (Schaef- 
fer, 1984, and Groth & Peebles, 1977): 



E^ a) En^ xi ' x ^ 

(a) 



(27) 



p-i 



where (a) is a particular tree-topology connecting the p 
points without making any loop, Qp™' is a parameter as- 
sociated with the order of the correlations and the topology 
involved, t a is a particular labeling of the topology, (a), 
and the product is made over the (p — I) links between the 
p points with two-body correlation functions. We show in 
Fig.f, taken from Valageas, Barber & Munshi (2003), the 
three topologies which appear within this framework for the 
5-point connected correlation function. 

Then, the minimal tree-model corresponds to the spe- 
cific case where the weights Q p a ^ are given by (Bernardeau 
& Schaeffer, 1992): 



Q ( P a) = n • 

vertices of (a) 



(28) 



where u q is a constant weight associated to a vertex of the 
tree topology with q outgoing lines. The advantage of this 
minimal tree-model is that it is well-suited to the compu- 
tation of the cumulant generating functions as denned in 
cq.(I6) for the density contrast Sr. Indeed, for an arbitrary 
real-space filter, -F(x), which defines the random variable s 



s= dxF(x)<S(x) and £ s = (s 2 ), 



(29) 



it is possible to obtain a simple implicit expression for the 
generating function, ip s {y) (see Bernardeau & Schaeffer, 
1992, and Jannink & Des Cloiseaux, 1987): 



ip s {y) = y j dx F(x 
r(x) = 



C[t(x)] 



r(x)C[r(x)] 



/dx'F(x')^^£[r(x')] 



(30) 
(31) 



where the function C„(t) is denned as the generating func- 
tion for the coefficients u v : 



(-1)" 



u„ t p with i/i = l. 



(32) 



Of course, the generating function, tp s (y), depends on the 
filter, F, which defines the variable, s. If the real-space fil- 
ter, -F(x), is close to a top-hat (this is actually the case for 
the smoothed density contrast, Sr, defined in eq.(14), where 
-F(x) is the top-hat of radius R normalised to unity), a sim- 
ple "mean field" approximation which provides very good 
results (Bernardeau & Schaeffer, 1992) is to integrate r(x) 
over the relevant volume V in eq.(31) with a weight F(x) 
and then to approximate r(x) by a constant r. This leads to 
the simple system (17)-(18) we have already encountered in 
Section 3.1. Therefore, in case the density field is described 
by such a minimal tree-model, the function f (t) introduced 
in Section 3.1 to parameterize the cumulant generating func- 
tion tp(y) can also be interpreted as a good approximation 
to the generating function C, v {t) of the vertices u p . 



3.4 Stellar model 

In a previous work (Valageas, Barber & Munshi, 2003) where 
we studied the properties of the shear, 7, we introduced an- 
other simple model for the many-body correlations of the 
density field which is well-suited to the computation of weak- 
lensing effects. This "stellar-model" is another case of the 
tree- models defined in (27), where we only keep the stellar 
diagrams (e.g., the graph (a) in Fig.l for the 5-point con- 
nected correlation) . Thus, the p— point connected correlation 
£ p of the density field can now be written as: 



£p(xi,..,x p ) = — y^]~[6(xi,xj). 

i=l j^i 



(33) 



The advantage of the stellar-model (33) is that it leads to 
very simple calculations in Fourier space. Indeed, eq.(33) 
reads in Fourier space: 

p 

(5(ki)..«5(k p )) c = ^ fo(ki k XII p (^)- ( 34 ) 

i — l j^i 

Of course, the Dirac factor, <5o(ki + ■■ + k p ), simply trans- 
lates the fact that the many-body correlations are invariant 
through translations. The coefficients S p and S p are related 
by: 



dt 3r 



Jf A 2 (fc)F 2 (fc.R) 



(35) 



where we introduced the Fourier transform, F(kR), of a 3-d 
top-hat of radius R: 



F(kR) 



f dx 
J V 



sin(kR) - (kR) cos(fcfl) 



However, in the following we shall use the simple approxi- 
mation: 



Sp. 



(37) 



p=i 



Alternatively, we may define the function ip(y) obtained 
from (17)-(18) and our choice of £(r) as the generating func- 
tion of the coefficients S p , rather than S p , through its Taylor 
expansion at y = 0. 
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3.5 Comments regarding our prescription versus 
the stable-clustering Ansatz and the halo 
model 

Here we must point out that our prescription where the 
p— point correlations are expressed from a tree-model but 
the coefficients S p (whence Q p a ^ in eq.(27)) are allowed to 
depend on scale and time is not a standard hierarchical 
model.* This is only a simple model which yields the angu- 
lar dependence of the many-body correlations, over a given 
scale, keeping the freedom to vary their amplitude with scale 
and time. Therefore, it is mainly a phenomenological tool 
which is not exactly self-consistent. Indeed, if we are given 
the two-point correlation and the coefficients S p at a given 
scale, applying the stellar model to all scales would only 
be fully self-consistent with constant S p . Otherwise, there is 
some ambiguity: given a configuration for the p— point corre- 
lation, which scale should be attached to this configuration 
so as to define the relevant parameters S p ? For our pur- 
poses, our prescription is sufficient and well-defined because 
smoothed weak-lensing observables like the smoothed con- 
vergence k s precisely select a "typical" length scale at each 
redshift along the line of sight: the radius of the smoothing 
cone V(z)6 s . In Fourier space, this is given by the typical 
wavenumber which we defined in eq.(22). Of course there re- 
mains a small ambiguity since we could as well have chosen 
1.5/(T>9 S ) but the associated uncertainty is smaller than the 
accuracy of our analytical model. 

The reason why we use such a two-step procedure is 
simply that we know the coefficients S p (such as the skew- 
ness S3) do depend on time and scale in CDM cosmologies. 
Therefore, it is necessary to keep such a freedom for the 
generating function <p(y), as described in Sect. 3.2. Then, as 
discussed below eq.(26), it happens that for some purposes 
the coefficients S p are not sufficient as one needs the an- 
gular dependence of the many-body correlations over some 
typical scale. This is why we introduced the two tree models 
presented in Sect. 3.3 and Sect. 3.4, which are understood 
to apply around some typical length scale as we explained 
above. 

As recalled in Sect. 3.2, our prescription in the highly 
non-linear regime for ip(y) is consistent with the stable- 
clustering Ansatz (Peebles 1980). However, some results 
from numerical simulations (e.g., Smith et al., 2002) sug- 
gest that the stable-clustering Ansatz is only approximate 
as continuing mergers yield deviations from this simple 
model. Such departures could be included within our ap- 
proach through the non- linear power-spectrum, P(k), and 
through an additional dependence on redshift and scale for 
the skewness, S3, in the non- linear regime. If required, one 
could simply choose another parameterization for the cu- 
mulant generating function <p(y) which would involve a spe- 
cific dependence on redshift and scale obtained from a fit 
to numerical simulations or from some other model for the 
density field. Indeed, our calculations do not depend on the 
parameterization (17)-(24) and we simply express the gener- 
ating function (y) of the convergence k s in terms of the 
generating function <p(y) of the density contrast. Therefore, 



* We would like to thank A. Taruya for calling our attention to 
the risks of confusion with more standard models. 



one simply needs to replace this function k s , z) by one's 
specific choice for the dependence on k 3 and z. 

In particular, an alternative to the stable-clustering 
Ansatz is provided by the "halo model" where the density 
field is described through a random distribution of dark mat- 
ter halos, modulated by the large-scale matter distribution 
(e.g., Seljak 2000, Peacock & Smith 2000). This allows one to 
include deviations from the stable-clustering Ansatz brought 
by the mergings and disruptions of these halos. However, 
as shown in Valageas (1999), note that the simplest halo 
model where all halos would have the same density pro- 
file (rescaled to their virial radius) is strongly inconsistent 
with numerical simulations as it yields £, p (R) ~ ft-wl 3 - 3 ) 
where /3 is the slope of the inner density profile of these 
halos (i.e. phaio(r) oc r _/3 ). Thus, the usual halo models in- 
volve a density profile which depends on the mass of the 
halo, through a concentration parameter. As noticed by 
Navarro et al. (1996) a good estimate is obtained by re- 
quiring the density within the core radius to scale as the 
density of the universe at the redshift when this mass scale 
turned non-linear. This can actually be seen as a way to 
include some features of the stable-clustering Ansatz within 
the halo model. Then, using such a dependence on M and 
n for the halo profiles (or a fit from simulations) one can 
match the observed non-linear power-spectrum (e.g., Smith 
et al. 2002) and possibly reproduce higher order correla- 
tions. Indeed, as noticed in Valageas (1999), if the low-mass 
tail of the halo mass function shows a power-law behaviour 
of the form n(M)dM/M oc M - 1 AM/M one recovers the 
stable-clustering Ansatz for 6 — > (in which case r/(M) also 
counts substructures within halos) so that one can probably 
obtain good results for higher order correlations by choosing 
a slope 9 <; 0. Then, one could compute from such a model 
the cumulants (<5^) c whence the generating function <p(y). 
Next, one can use all methods presented in this paper with 
this new function <p(y). 

However, it is not clear whether such a halo model can 
be made fully self-consistent since higher order correlations 
may be increasingly influenced by the substructures present 
within halos (see also Valageas 1999). Moreover, this model 
cannot reproduce the behaviour of the density field over low- 
density regions like voids and filaments (i.e. outside virial- 
ized halos). This prevents the computation of the full pdfs 
V(S R ) and V(k) for S R <> and k <> 0. 

Therefore, in this work we shall only consider the sim- 
ple parameterization (19)-(24) described in Sect. 3.2. Indeed, 
detailed comparisons with numerical simulations, presented 
below in Sect. 6, show that it already provides good pre- 
dictions for the pdf V(k) of the convergence. Moreover, this 
model has the advantage of a great simplicity and it auto- 
matically shows the right behaviour on quasi-linear scales. 



4 COMPUTATION OF THE PDFS 

We now compute the pdf V(k) of the convergence following 
the method developed in Valageas (2000a, b). From eq.(7) or 
eq.(9) we express the cumulants of the smoothed normalised 
convergence (&J ) c in terms of the many-body correlations of 
the density contrast. Next, after resummation of this se- 
ries of cumulants we obtain the generating function ip ks of 
the convergence as in eq.(16) which yields the pdf V{k s ). 
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Since this method has already been used in previous works 
for the convergence (Valageas 2000a, b), the aperture- mass 
(Bernardeau & Valageas 2000) and the shear (Valageas, Bar- 
ber & Munshi 2003) we only briefly recall the main steps of 
this derivation, in order to show where the various approxi- 
mations one can introduce within this framework come into 
play. 

We first present in Sect. 4.1 and Sect. 4.2 two simple 
approximations which allow us to compute V[n) without 
any assumption about the properties of the density field. 
They provide two simple expressions for (pa s (y) in terms 
of <p(y) which can be used with our parameterization de- 
scribed in Sect. 3.2 or with any alternative prescription for 
<p{y) or V(5r). This allows us to consider in Sect. 4.3 a log- 
normal model for V{8r). Finally, in Sect. 4.4 and Sect. 4.5 
we present exact calculations within the framework of both 
the minimal tree-model and the stellar-model introduced in 
Sect. 3.3 and Sect. 3.4, in order to check in two explicit 
cases the model-independent approximations described in 
Sect. 4.1 and Sect. 4.2. 

4.1 Spherical-cell approximation 

Following Valageas (2000a, b), see also Bernardeau & Valageas 
(2000), we obtain from eq.(7) for the cumulant of order p of 
the smoothed normalised convergence k 3 : 

JO J-oc i=2 Jo i=1 na ° 



X £p 







X2 



Xp 



(38) 



Here we used the fact that the correlation length (be- 
yond which the many-body correlations are negligible) is 
much smaller than the Hubble scale, c/H(z) (where H(z) 
is the Hubble constant at redshift z). Although the points 
(Xi,C^i) cover a cylinder of radius V9 S and length L (with 
L — > oo in eq.(38)) rather than a sphere, we may approxi- 
mate the integral over the p— point correlation £ p as: 

n^/'nS* 6 '-^" 1 ' (39) 



- oo - o J 



in a fashion similar to £ p = S p £ 2 > see eqs.(16), (26). Here 
the coefficients S p (z) are evaluated at the wavenumber k s 
defined in eq.(22), associated with the radius T>{z)9 s of the 
cylinder at redshift z, and we defined: 

/■oo p6. 

I, 



lK{t)= L dx L 



This quantity can also be written as (Valageas 2000b): 
dfc A 2 (k,z) 



k 



w(kve s y 



(41) 



where the filter W was defined in eq.(10). Thus, using eq.(39) 
we can write the cumulants (38) as: 



dx w p S p II 



p-i 



(42) 



Then, using the expression (16) for the cumulant generating 
function tp ks (y) associated with the smoothed normalised 
convergence k s we obtain: 



(43) 



where we used the resummation (16) for the coefficients S p 
and we introduced the variance: 



U s = (£) = / 

Jo 



In- 



(44) 



Of course, eq.(44) is exact, within the small-angle approxi- 
mation, and the pdf V{k s ) is finally obtained from the in- 
verse Laplace transform (15): 



V(Ks) 



s: 



dy 



(45) 



We shall refer to the approximation (39)-(43) as the "spherical- 
cell approximation" because it is based on the use for the 
average of the many-body correlations over a cylinder (39) 
of their average over a spherical cell: 



(Cp)cyl. ^ (Cp)sphcr. _ g 
spher 



(46) 



In this equation only the notation (..) is a spatial aver- 
age, over spherical cells or cylinders. This approximation 
should not be confused with an approximation based on 

some "spherical dynamics". It merely assumes that the de- 

— — p— i 

pendence on geometry of the ratios Cp/^2 can be ne- 
glected. The advantage of this approximation is that it can 
be applied to any model for the density field. Indeed, it 
does not involve the detailed behaviour of the many-body 
correlations £ p (xi, .., x p ): we only need the pdf V(8r) (or 
the associated generating function <p{y)) of the density con- 
trast within spherical cells. In particular, we shall see in the 
following sections that eq.(43) can be recovered from more 
specific models. As recalled in Sect. 3.3, the drawback of this 
approximation is that it cannot be extended to compensated 
filters like those used for the aperture-mass or the shear. 



4.2 Mean-redshift approximation 

We see in eq.(43) that the projection of the 3-d density field 
onto the 2-d convergence is described by a simple integra- 
tion along the line of sight of the generating function <p(y) 
(within the framework of the spherical-cell approximation). 
We wrote explicitly in eq.(43) the dependence on redshift 
of <p(y) along the line of sight, which follows the evolution 
of the wavenumber probed at redshift z, see eq.(22), as well 
as the growth of density fluctuations from the linear to non- 
linear regime, see eq.(24). As noticed in Valageas (2000a, b), 
we may approximate this integral along the line of sight by a 
mean value. Since both Xs£,k. s /Ik and wI K /£,k B are typically 
close to unity one may simply use: 



(47) 



which still obeys the constraint ipk s (y) = —y 2 /2 + ... at y — 
0, see eq.(16). In eq.(47), we take for the generating function 
<p{y) its value at the typical redshift z* (and wavenumber 
k s (z*)) which we define as the location of the maximum 
of the selection function w(z,z s ). This yields z* ~ z s /2. 
We shall refer to this approximation (47) as the "mean- 
redshift approximation" because it replaces the projection 
onto two dimensions by a typical value along the line of sight. 
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Eq.(47) actually means that the pdf V(k 3 ) of the smoothed 
normalised convergence is directly given by the pdf V(5r) of 
the density contrast at the typical scale l/k s (z*) and redshift 
z* probed by the observation: 



V(k s ) ~ V{S R -> k s ,£ 2 -> £ ft , ;&„(;?»), 2„). 



(48) 



Eq.(48) clearly shows that the statistical properties of the 
convergence provide an efficient probe of the density field. 
In particular, we can hope to measure the departures from 
Gaussianity brought by the non-linear gravitational dynam- 
ics from V{k s ). 



4.3 Log-normal approximation 

It is clear that both the spherical-cell approximation and 
the mean-redshift approximation presented in Sect. 4.1 and 
Sect. 4.2 are model-independent. Indeed, they do not assume 
any specific behaviour for the generating function ip(y) or 
the many-body correlations £ p . However, the accuracy of 
these simple approximations may depend on the properties 
of the density field. For instance, if tp{y) or the selection 
function w strongly evolve with redshift the average used 
for ip(y) within the mean-redshift approximation (47) may 
be too inaccurate. Nevertheless, within the framework of 
the mean-redshift approximation (47)- (48) one may directly 
use any model for the pdf of the density contrast to estimate 
V(k s ), with no need to compute the generating function <p(y) 
itself. Thus, Taruya et al. (2002) used this mean-redshift 
approximation to compute the pdf V{k s ) from a log-normal 
distribution for the density contrast. That is, they used for 
V(8r) in eq.(48) the expression: 



T\n{8 R ) 



(l + ^) v / 27rln(l + C 2 



x cxp — 



ln 2 [(l + <Wi+£ 2 
21n(l + £ 2 ) 



(49) 



As noticed by Taruya et al. (2002) the log-normal pdf is a 
simple example which violates the stable-clustering Ansatz. 
However, as we explained in Sect. 3.5 our method does not 
need the stable-clustering Ansatz to be valid so that the 
log-normal pdf is fully consistent with the approach recalled 
in Sect. 4.1 and with eqs.(47)-(48). One advantage of this 
model over the use of the generating function ip(y) obtained 
from cqs.(17)-(18) is that one does not need to compute 
the inverse Laplace transform (15) and there is one fewer 
parameter, indeed the log-normal pdf (49) only involves the 
variance £ 2 . On the other hand, we can expect the additional 
parameter ft (or S3) which enters our parameterization (19) 
and allows us to take into account the dependence of the 
skewness on the slope of the power-spectrum to improve the 
accuracy of this prescription over the log-normal approxi- 
mation. 

For instance, as noticed in Bernardeau (1994), in the 
quasi-linear regime the log-normal approximation is quite 
good for n = — 1 but worsens for different power-spectra. We 
may note here that in the quasi-linear limit our prescription 
(17)-(19) actually yields back the log-normal pdf for S3 — > 3, 
where f (r) — > e~ T — 1 (e.g., App.A of Valageas, 2002). From 
eq.(21) this corresponds to n = —1.14 which explains why 
the log-normal pdf agrees with numerical simulations for 



n ~ — 1 in the quasi-linear regime. However, we can expect 
significant discrepancies in other cases where S3 is much 
larger than 3 (e.g., for n ^ — 2 or in the highly non- linear 
regime) . 



4.4 Minimal tree-model approximation 

In the spherical-cell approximation presented in Sect. 4.1 
we estimated the average over cylinders of the correlations 
£ p by the simple approximation (39). As discussed above, 
the advantage of this approach is that it does not require 
much information about the behaviour of the density field. 
However, one may wonder what is the actual accuracy of this 
simple estimate. To tackle this point we need to evaluate 
exactly the l.h.s. in eq.(39) for some specific models and 
compare the result with the r.h.s. approximation. We could 
also compute both quantities directly in N-body simulations. 
Note that the pdf P(k s ) measured in numerical simulations 
does not allow us to test the approximation (39) itself since 
it also involves the parameterization used for the generating 
function tp(y). Therefore, it is interesting to compute the 
pdf V{k s ) in a more rigorous way, without using eq.(39), for 
some specific cases. 

Obviously, to do so, we need the detailed behaviour of 
the p— point correlation £ p (xi, ..,x p ). One simple example 
where we can sum up the cumulants (k?) c so as to derive 
the generating function ip fia (y) is the minimal tree-model 
recalled in Sect. 3.3. This case has already been studied in 
Valageas (2000b) and Bernardeau & Valageas (2000). Here 
we shall simply recall how we can recover the spherical-cell 
approximation (39) within this framework. As seen from 
eq.(38), in order to make some progress we need to eval- 
uate the quantities: 



(50) 



Then, as noticed in Valageas (2000b), if the 3-d correlations 
£ p obey a tree-model as in eq.(27) the 2-d correlations lo p 
exhibit the same tree-structure: 

cu p (S 1 ,..,S p ;z)=J2Q i p a) J2Il UJ2 ^'^'^ ( 51 ) 

(ct) t a p-1 



with: 

W2(l?l,1?25«) = 7T 



dfe A 2 (fc,z) 



k 



k 



Jo(fcD|i?i-0 2 |), (52) 



where Jo is the Bessel function of order 0. Next, in the case of 
a minimal tree-model (28) we can perform the resummation 
(30)-(31) for the 2-d correlations lo p , since the latter obey 
the same minimal tree-model from eq.(51). This yields (see 
Bernardeau & Valageas 2000 for details): 



Vft, (y) 



d X -j- Vcyl. 



. Ik 

yw—\z 



(53) 



where we introduced the 2-d generating function ip cy i. associ- 
ated with the 2-d correlations ui p , given by the resummation: 



fcy\.{y) 



Ad 



C[t(*)] 



t(0)£[t(0)] 



(54) 
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T(0) 



f 6s d#' 



(55) 



Here we introduced the angular average LO2 of the 2-d corre- 
lation lo 2 , which is actually equal to the quantity I K defined 
in eqs.(40)-(41): 



u>2(z) 



-r 

Jo 



^2 ^2 w 2 (i?i,i?2;«) = /«(«), 



(56) 



and we substituted the notation I K to ZJ2 into eq.(53). Note 
the similarity of eq.(53) with eq.(43) obtained from the 
spherical approximation. Moreover, we see that for a min- 
imal tree-model the projection from three dimensions onto 
two dimensions does not only yield the integration along the 
line of sight apparent in both eq.(53) and eq.(43). It also en- 
tails the change from ip(y), associated with 3-d spherical 
cells, to (p C yi. , associated with a 2-d top-hat (after we inte- 
grated along the longitudinal direction, see eq.(50)). Thus, 
this projection effect was actually neglected by the spherical- 
cell approximation (39). However, if we use the mean-field 
approximation discussed below eq.(32) (i.e. r(x) is approx- 
imated by a constant over the relevant volume of integra- 
tion), both for 3-d spherical cells (which amounts to the 
approximation £(t) ~ C^( r )) and for the 2-d top-hat of an- 
gular radius 9 S which appears in eqs.(54)-(55), we obtain 
fcy\.{y) — f(y) and we recover eq. (43). Therefore, within the 
framework of a minimal tree-model for the many-body cor- 
relations the spherical approximation presented in Sect. 4.1 
can be interpreted as the usual mean-field approximation. 
This is not surprising since the latter approximation actu- 
ally neglects the geometry of the filter -F(x) associated with 
any random variable s, see eqs.(29)-(31), so that the coeffi- 
cients S p must be the same for spherical cells and cylinders. 

The mean-field approximation has already been shown 
to provide very good results in the case of the 3-d top-hat 
(Bernardeau & Schaeffer 1992) hence we can expect a sim- 
ilar accuracy for the 2-d top-hat. This suggests that the 
spherical-cell approximation (39) should be quite accurate 
for most models of the density field, even beyond the class 
of minimal tree-models. 



4.5 Stellar-model approximation 

Finally, we consider the second explicit model presented in 
Sect. 3: the stellar model introduced in Sect. 3.4. The cal- 
culation of the generating function (pa s (y) can again be per- 
formed explicitly and we can also recover the result (43) 
obtained from the spherical-cell approximation through an- 
other mean field approximation. Moreover, we shall compare 
the exact results obtained from this approach with those of 
the spherical-cell approximation and we shall check in the 
Figures that both predictions are indeed very close, in agree- 
ment with the discussions above. The stellar model (33)- 
(34) was already introduced in Valageas, Barber & Munshi 
(2003) to compute the pdf of the smoothed shear compo- 
nents, 71s and 72s, as well as the pdf of the smoothed mod- 
ulus, 7 S = |7is + 172s I, hence we shall only recall here the 
main steps of the calculation since it proceeds in the same 
fashion for the smoothed convergence, n s , which we study in 
this paper. Within this framework it is convenient to work 



in Fourier space, hence we write the cumulant («?) c as: 




n** / n 



dkj W{k±jVe s ) 



S P fo(ki + + k p )P(k 2 )..P(k P ), (57) 



where we used eq.(9) and eq.(34), and the fact that the 
correlation length is much smaller than the Hubble scale. 
Thus, eq.(57) is equivalent to eq.(38) which was written in 
real space, where we used the stellar model (34) for the 
correlation (<S(ki)..<S(k p )} c . Next, using the usual small-angle 
approximation (i.e. P(kj) ~ P(k±j)), we can perform the 
integration over \2, ■■, Xp and fc||i, ••, k\\ p , which yields: 



x S p fo(kxi + .. + k ±p )P(k ±2 )..P(k ±p ). (58) 



Then, using the standard exponential representation of the 
Dirac distribution (see Valageas, Barber & Munshi 2003), 
integrating over the angles of the transverse wavenumbers 
kxi, .., kx p and over |kj_i|, we obtain: 



(k p s)c= I** d X w p f dt2tS p P K ~\ 
Jo Jo 



where we introduced the quantity 
dfc A 2 (k,z 



I K *(t, z) = TV 



W(kV6 s ) Jo(tkV6 s 



(59) 



(60) 



We used an asterix "*" in the notation I K *(t, z) in order to 
distinguish this quantity from I K [z) introduced in eqs.(40)- 
(41) and to recall that 7 K * appears within the stellar model. 
Finally, we can resum the cumulants (kJ) c , using the re- 
summation (16) where we replace S p by S p (as discussed in 
Sect. 3.4 this may also be taken as the definition of ip(y)), 
and we get: 



(61) 



<Pk.{v)= I* d X f &t2t^ip(yw 1 -^;. 
Jo Jo 1k * \ 

Thus, we again obtain a result very similar to the spherical- 
cell approximation (43). As for the minimal tree-model, we 
see through this second exact calculation that the projec- 
tion effects also lead to a modification of the 2-d generat- 
ing function associated with a given redshift z as compared 
with <fi(y). For the minimal tree-model, this led to the 2-d 
generating function ip C yiXv) obtained in eqs.(54)-(55). For 
the stellar model, this projection effect is described by the 
integration over t of the generating function ip which ap- 
pears in the r.h.s. in eq.(61). Next, we can again recover the 
spherical-cell approximation by using a suitable average for 
the integration over t. Indeed, noticing that J dt 2t = 1 

and J dt 2tI K * = I K , we may use the approximation: 



f dt 2t if ( 
Jo Ik * \ 



yw — ;z 



In 



yw — ;z 



(62) 



which yields back the spherical-cell result (43). We shall 
check through numerical computations that the approxima- 
tion (62) is indeed very good, which gives a further justifi- 
cation for the simple spherical-cell approximation (43). 
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5 CONVERGENCE STATISTICS FROM 
NUMERICAL SIMULATIONS 

Our results for the pdfs of the convergence and the higher- 
order moments established from the above analytical models 
are compared in this paper with our results obtained from 
numerical simulations. The numerical method we have used 
was originally devised by Couchman et al. (1999) who devel- 
oped a code for computing the 3-dimensional shear matrices 
at locations within the simulation volumes. We have ap- 
plied this code to the simulations of the Hydra Consortium^ 
produced using the 'Hydra' JV-body hydrodynamics code 
(Couchman, Thomas & Pearce, 1995). 

In this pilot study, in which analytical predictions are 
made and compared with the results from numerical simu- 
lations for arbitrary cosmological scenarios, we have chosen 
readily available simulation data. Details of the two cos- 
mological simulations used, LCDM and OCDM, are sum- 
marised in Table 1, in which the angular size of the survey, 
^survey i refers to the total angular size of the field of view 
and corresponds to completely filling the front face of the 
redshift 1 simulation volume. For source redshifts greater 
than 1, the periodicity of the particle distributions was used 
to allow lines of sight beyond the confines of the simulation 
volumes to be included. This is unlikely to adversely affect 
the statistics for the angular scales of interest, i.e., l'.O to 
8'.0. In addition, there is, of course, no affect for the source 
redshifts less than or equal to 1. 

Both cosmologies contained 86 3 dark matter particles of 
mass 1.29 x lO 11 /). -1 solar masses each, where h expresses the 
value of the Hubble parameter in units of 100 km s _1 Mpc -1 . 
A variable particle softening, whose value reflected the den- 
sity environment of each particle, was used. The minimum 
value (in box units) of the softening (for particles in the 
densest environments) was 0.0007(1 + z), where z is the red- 
shift of the particular simulation volume. Each simulation 
volume had comoving side-dimensions of 100/t -1 Mpc and 
to avoid obvious structure correlations, each was arbitrarily 
translated, rotated and reflected about each coordinate axis. 

The full procedure for specifying the coordinates for 
the lines-of-sight, the locations within the simulations for 
the computations of the shear and the procedure for the 
combination of the 3-d matrices to obtain the final Jaco- 
bian matrices is described in detail by Barber (2002) . In the 
present work, a total of 455 x 455 lines of sight were used 
and 300 regularly-spaced evaluation locations were specified 
along each line of sight in each simulation volume. In the 
LCDM cosmology, the full field of view was 2°. 6 x 2°. 6 and 
in the OCDM cosmology, 2° .8 x 2° .8. The angular resolu- 
tion in the LCDM cosmology was 0'.34, which equates to 
the minimum value of the particle softening at the optimum 
redshift, z = 0.36, for lensing of sources at a redshift of 1. 
In the case of the OCDM cosmology, the angular resolution 
was 0'.37. 

A total of 14 source redshifts were selected in each cos- 
mology to give good statistical coverage of the redshifts of 
interest (see Table 2). Each complete pass through all the 
simulation volumes and for each source redshift was per- 
formed a total of N — 10 times and the values for the con- 
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Table 1. The parameters used in the two cosmological simula- 
tions, r is the power spectrum shape parameter, f! m is the matter 
density parameter, S7a is the vacuum energy density parameter, 
as is the normalisation on scales of 8/i _1 Mpc, 8 rC s is the angular 
resolution and 9 2 UTVey is the angular size of the complete field 
throughout the simulations. 
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Figure 2. Variance of the smoothed convergence field k„ (k 2 ) 1 / 2 , 
as a function of the source redshift z s . The smoothing angle 8 S is 
fixed at 1 arcminute in left panel and at 4 arcminute in the right 
panel. Lines correspond to the analytical prediction (44) and data 
points represent results from numerical simulations (filled circles 
for LCDM and open circles for OCDM). Error bars are computed 
from scatter among various realizations. 



vergence were smoothed on the different angular scales using 
a top-hat filter. The computed values for the pdfs and the 
higher-order moments from each of the runs in each cosmol- 
ogy were averaged so that the errors on the means of la/^/N 
for each statistic were determined. 



6 RESULTS 

Comparison of our analytical results against numerical simu- 
lations can be divided in two different categories. In addition 
to comparing the analytical results for lower order moments 
we also compute the whole pdf V(k 3 ) of the smoothed con- 
vergence k s for a wide range of smoothing angles and source 
redshifts. This allows us to assess the accuracy of the simple 
parameterization described in Sect. 3.2 and of the various 
methods presented in Sect. 4. 



6.1 Amplitude of the convergence 

We first study here the amplitude of the smoothed conver- 
gence as a function of angle and source redshift, for both 
LCDM and OCDM cosmologies. Thus, we show in Fig. 2 
the dependence on the source redshift z 3 of the rms con- 
vergence (kJ) 1 ' 2 for two smoothing angles, a = 1' and 
9 s — 4' . The lines show the analytical prediction (44) while 
the data points are the results from numerical simulations. 
The rms convergence increases with the source redshift, z s , 
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Table 2. The redshifts of the sources in each cosmology were selected to give good statistical coverage throughout the overall range. 
In each cosmology, they corresponded to the simulation box redshifts and were chosen to be close to redshifts of 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 3.0 and 3.5. In this paper we refer to the source redshifts as these approximate values, although in the 
determination of our results, the actual redshift values, listed below, were used. 
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Figure 3. Variance of the smoothed convergence field re s , (re 2 ) 1 / 2 , 
as a function of smoothing angle 8 S for three source redshifts z s . 
These are z B = 0.5 (lower curve), z s = 1 (middle curve) and 
z s = 1.5 (upper curve). Errorbars are computed from scatter 
among various realizations. Different panels correspond to differ- 
ent cosmologies as indicated. 
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Figure 4. Minimum value |re m i n | of convergence re as a function 
of source redshift z s . Note that re m ; n is independent of smoothing 
angle 9 B . 

according to the (cosmology-dependent) rate of formation 
of structure and the location of massive structures for their 
gravitational lensing effects, as described by Barber et al., 
2000. This makes the pdf broader. The variance is smaller 
for the OCDM case mainly because the normalization, as, 
of the power-spectrum, P(k), is smaller (see Table 1). In a 
similar fashion, (k 2 } 1//2 decreases for larger smoothing an- 
gles which probe larger scales, see eq.(22), where the am- 
plitude of the density fluctuations is smaller. We obtain a 
reasonable agreement with the results from numerical sim- 
ulations which shows that the prescription from Peacock & 
Dodds (1996) is sufficient to reproduce the non- linear power- 



spectrum over the regime probed by these angles and source 
redshifts. Note indeed that eq.(44) for the variance £„ s only 
involves the Born and small-angle approximations, as well 
as our prescription for P(k) (here taken from Peacock & 
Dodds, 1996). 

Next, we plot in Fig. 3 the dependence on the smooth- 
ing angle, 8 a , of the rms convergence, (k; 2 } 1 / 2 , for the three 
source redshifts, z s — 0.5, 1 and 1.5. The redshift, z s , cor- 
responding to the curves increases from bottom to top, in 
agreement with Fig. 2. Of course, as in Fig. 2 we can check 
that the variance grows at smaller angular scales which 
probe deeper within the non-linear regime of gravitational 
clustering. Consistently with Fig. 2 we again obtain a reason- 
able agreement with the results from numerical simulations, 
although we seem to overestimate somewhat the variance for 
the OCDM cosmology. However, this overestimate is within 
the errors associated with iV— body simulations and the Pea- 
cock & Dodds (1996) fit to the non- linear power-spectrum. 
We note that the Peacock & Dodds (1996) fit underestimates 
the more recent simulation results of Smith et al. (2002) in 
the range 1 < A 2 (fc) < 50 for the LCDM and OCDM cos- 
mologies, whilst it overestimates their numerical data for 
higher wavenumbers. Our overestimate might be due to the 
fact that the slope of the power, A 2 (fc), is somewhat steeper 
for OCDM (see Smith et al., 2002, Figure 15), which means 
that the relative contribution from high wavenumbers should 
be larger; whence the Peacock & Dodds (1996) fit would be 
expected to overestimate the simulation results by a larger 
amount for the OCDM cosmology than the LCDM. 

Finally, we show in Fig. 4 the lower-bound K m i n of the 
convergence k as a function of the source redshift z s , see 
eq.(5). Of course, |ft m in| increases at higher z a with the 
length of the line of sight. It is interesting to compare |fv m i n | 
shown in Fig. 4 with the rms convergence (kj) 1 '' 2 shown in 
Fig. 2. Indeed, it partly describes the deviation of the pdf 
V(k s ) from the Gaussian. For |ft m m| ^> ( K s} 1//2 the lower- 
bound Kmin has a strong influence on the shape of the pdf 
which has to be significantly different from the Gaussian. 
On the contrary, for |ft m in| ^> (k 2 ) 1 ^ 2 this lower-bound is 
located very far in the low convergence tail of the pdf which 
can therefore look roughly similar to a Gaussian (the asym- 
metry is weaker). Comparing Fig. 4 with Fig. 2 we see that 
the pdf will be more asymmetric for low source redshifts z a . 
This trend also follows from the fact that low redshifts also 
probe the late stages of gravitational clustering where the 
density field has evolved farther from the Gaussian initial 
conditions. 
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Figure 5. The typical comoving wavenumber k s (in units of 
MVlpc -1 ) probed by the smoothed convergence K a as a function 
of the redshift z along the line of sight, from eq.(22). The up- 
per pair of curves corresponds to the smoothing angle 8 S = 1' 
whereas the lower set of curves corresponds to the smoothing an- 
gle B =4'. Note that the redshift z is not the source redshift 
z a but the redshift along the line of sight (0 < z < z a ) of the 
intermediate lensing structures. 
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Figure 6. The power per logarithmic wavenumber interval 
A 2 (k s (z), z) as a function of redshift z along the line of sight. The 
upper pair of curves corresponds to the smoothing angle 6 S = 1' 
while the lower set of curves corresponds to the smoothing angle 
6> s = 4'. 



6.2 The regime of gravitational clustering probed 
by weak-lensing 

As the smoothing angle 8 3 and the source redshift z s vary the 
convergence n s probes different scales and different regimes 
of gravitational clustering. This could allow one to derive 
some information about the physics of the gravitational dy- 
namics in the expanding universe from future weak-lensing 
surveys, in addition to the measure of the main cosmological 
parameters. 

First, because of the dependence on redshift of the angu- 
lar distance T>(z), eq.(4), the typical comoving wavenumber 
k 3 (z) probed by the smoothed convergence k s varies with 
the redshift z along the line of sight, see eq.(22). Thus, we 
display in Fig. 5 this typical comoving wavenumber k s as a 
function of z, for both LCDM and OCDM cosmologies and 
both smoothing angles 6 S = l' and 9 a = 4'. Note that the 
redshift z is not the source redshift z s but the redshift along 
the line of sight (0 < z < z s ) of the intermediate lensing 
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Figure 7. Spectral index n(k s (z)) as a function of redshift z 
along the line of sight. The upper pair of curves corresponds to the 
smoothing angle 8 a = 1' while the lower set of curves corresponds 
to the smoothing angle 8 S =4'. 

structures which give rise to the deflection of the light rays. 
As seen from eq.(22), the wavenumber k s is actually propor- 
tional to 1/9 S as larger smoothing angles probe larger scales 
and smaller wavenumbers. We see that for 6 a = 1' the typ- 
ical comoving wavenumbers probed by the convergence are 
of order 1 to 10/iMpc -1 which corresponds to scales 0.1 to 
l/i _1 Mpc. These are the scales of present galaxies and clus- 
ters. In particular, this means that for such angular scales 
weak lensing mainly probe the intermediate regime of grav- 
itational clustering. 

This is clearly apparent in Fig. 6 which shows that 
the power per logarithmic wavenumber interval A 2 (k s (z), z) 
at the scales probed by weak-lensing typically runs from 
A 2 (k s (z),z) ~ 1 up to A 2 (k s (z),z) ~ 400. Note that the 
properties of the density field show a fast evolution in this 
transition regime. This entails a critical test of our simple 
parameterization described in Sect. 3.2, which must be able 
to follow the evolution of gravitational clustering from lin- 
ear to highly non-linear scales. Moreover, we can suspect 
the mean-redshift approximation (47) to fail to reproduce 
the results from numerical simulations with a high accuracy 
in this transition range. Indeed, the value (p(y;z*) of the 
generating function at the typical redshift z* may be a low- 
accuracy approximation to the mean (43) of the generating 
functions (p(y; z) which characterize the density fluctuations 
encountered along the line of sight, as the latter show a sig- 
nificant evolution from z = up to z — z s . This point will 
appear clearly in Sect. 6.3 where we discuss the skewness of 
the convergence shown in Fig. 9. 

Finally, we show in Fig. 7 the slope n of the linear 
power-spectrum at the scales probed by the smoothed con- 
vergence along the line of sight. We see that it runs from 
n ~ —2.5 (for small angles and redshifts) up to n ~ — 1 
(for large angles and redshifts). Together with Fig. 6, this 
shows that gravitational weak-lensing actually probes a wide 
range of physical conditions for the underlying density field. 
This requires the use of flexible models which can cover any 
hierarchical power-spectrum and follow the gravitational dy- 
namics from the quasi-linear regime up to the highly non- 
linear regime, for arbitrary cosmological parameters. Fortu- 
nately, we shall see in the following sections that the sim- 
ple model presented in Sect. 3.2, which we use in this ar- 
ticle (except for the log-normal approximation introduced 



© 0000 RAS, MNRAS 000, 000-000 



14 Barber A. J. et al. 



LCDM 

OCDM 




Figure 8. Skcwness S 3 3 of the smoothed convergence field re 



a function of the source redshift z s . The left panel corresponds 
to the smoothing angle S = 1' whereas the right panel is for 
9 a =4'. The lines show our analytical prediction (63) from the 
spherical-cell approximation while data points represent results 
from numerical simulations. Error bars are computed from the 
scatter among various realizations. 




Figure 9. Skewness S3 3 of the smoothed convergence field 
a function of smoothing angle 9 S (in arcminute) for three source 
redshifts z s . These are z s = 0.5 (upper curve), z s = 1 (middle 
curve) and z s = 1.5 (lower curve). Different panels correspond 
to different cosmologies as indicated. The solid lines show the 
spherical-cell approximation (63) whereas the dotted lines corre- 
spond to the mean-redshift approximation (65). The predictions 
from the stellar model are almost identical to spherical-cell ap- 
proximation and are not reproduced here for clarity. 



in Sect. 4.3), is able to recover the results obtained from 
numerical simulations with a good accuracy over all ranges 
of interest. On the other hand, from an observational point 
of view, we note that this feature makes it somewhat more 
difficult to derive precise constraints on the non-linear grav- 
itational dynamics from weak-lensing surveys. Indeed, it is 
difficult to extract one precise scale and gravitational regime 
from the observed smoothed convergence k 3 . 

6.3 The skewness S3 3 of the smoothed 
convergence 

The simplest measure of the progress of gravitational clus- 
tering from small Gaussian initial conditions is provided by 
the third-order moment of the density field (<5fj). As dis- 
cussed in Sect. 3.1, it is actually more convenient to study 
the skewness S3 defined by the ratio S3 = (5fj)/(<5fj) 2 , see 
eq.(16). Since the convergence is linear over the density field, 
see eq.(l), it is natural to consider also the skewness S3 3 



of the smoothed convergence re s , defined in the same fash- 
ion. Thus, we show in Fig. 8 the dependence of S3 3 on the 
source redshift z s . For clarity, we only plot the result from 
the spherical-cell approximation, which yields from eq.(42): 



spherical-cell : £3 * = 



(ftj) / d X W 3 S3{Z)H 



[Jdx 



w 2 I K 



and: 



S? s = 



ML 



(63) 



(64) 



We consider both LCDM and OCDM cosmologies, and the 
two smoothing angles 8 S = 1' and 8 S — 4' . We can check in 
the figure that the skewness S3 3 decreases for larger source 
redshift z s . This is partly due to the fact that higher red- 
shifts probe earlier stages of gravitational clustering where 
the density field is closer to Gaussian and its skewness 53 
is also smaller. However, most of this decrease is due to the 
sum along the line of sight over successive mass sheets. In- 
deed, the total convergence re is the integral over the whole 
line of sight of the random contributions associated to the 
successive lens planes at redshifts < z < z s , see eq.(l). 
Then, this sum of random variables tends to make the total 
signal closer to Gaussian through the central limit theorem. 
This can also be seen as follows. Within the mean-redshift 
approximation (47) we have for the skewness S3 3 of the nor- 
malized smoothed convergence re s and for S3 3 : 



mean-redshift 



S 3 (z* 



and S? 



S3(z*) 



(65) 



Thus, we see from eq.(65) and from the increase of |re m i n | 
with the source redshift z s shown in Fig. 4, that even for 
a constant skewness £3(2) for the density field we would 
obtain a decrease with the source redshift for the skewness 
S3 3 of the smoothed convergence (while the skewness S3 3 of 
the normalized smoothed convergence re s would remain con- 
stant within the mean-redshift approximation). Of course, 
as noticed above, this trend is actually reinforced by the 
fact that higher redshifts also probe density fields which are 
closer to Gaussian (i.e. Ss(z) declines at larger z). From the 
behaviour of the skewness S3 3 we can already infer that the 
pdf V(Ks) of the convergence will look closer to Gaussian 
for higher source redshifts. 

We can see in Fig. 8 that our analytical prediction (63) 
shows a good agreement with the results from numerical 
simulations. This means that our simple parameterization 
(21)-(24) for the skewness of the density field provides a 
reasonable estimate over the range probed by weak-lensing. 
Indeed, the prediction obtained from the stellar model ap- 
proximation is almost identical to the result computed from 
eq.(63) which means that the spherical-cell approximation 
is very good (at least for low-order moments like the skew- 
ness) and that the comparison with numerical simulations 
in Fig. 8 mainly tests our simple parameterization (21)-(24) 
for the skewness of the density field. Note the rather small 
dependence on S7a of the skewness S3 3 . On the other hand, 
the skewness is well-known to show a strong dependence on 
the cosmological parameter Sl m , which explicitly appears in 
re m i n , see eq.(5). This provides a good tool to measure n m 
from weak-lensing surveys, see Bernardeau et al. (1997). 

Next, we plot in Fig. 9 the dependence on the smoothing 
angle 9 3 of the skewness S3 3 for the three source redshifts 
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z a — 0.5, 1 and 1.5. In agreement with Fig. 8 the redshift z 3 
corresponding to the curves increases from top down to bot- 
tom. The solid lines are the spherical-cell approximation (63) 
while the dotted lines are the mean-redshift approximation 
(65). We can see in Fig. 9 that the skewness S^ s decreases 
for large smoothing angles 6 S . This is due to the fact that 
large angles probe large scales which are closer to the lin- 
ear regime, whence their skewness S3 is smaller. This trend 
is also reinforced by the fact that for smaller wavenumbers 
k the local slope n of the linear power-spectrum increases 
(for CDM-like power-spectra like those we study here) which 
leads to a smaller skewness 5*3 in both the linear and highly 
non- linear regimes, see eq.(21). We can check that our pre- 
diction (63) agrees reasonably well with the numerical simu- 
lations which shows that our simple parameterization (21)- 
(24) works fairly well. The results obtained from the stellar 
model approximation are almost identical to the solid lines 
which means that the spherical-cell approximation is quite 
good. In fact, the inaccuracy is dominated by far by the 
parameterization (21)-(24) rather than by the spherical-cell 
approximation. 

On the other hand, we note that the mean-redshift ap- 
proximation (65) shows some discrepancies with the simula- 
tion results. In particular, it yields a characteristic step-like 
profile with a sharp decrease for larger angles. This step cor- 
responds to the transition to the highly non-linear regime 
where the skewness suddenly shows a sharp variation, as 
it evolves from S^ L up to S^. The plateau at small an- 
gles corresponds to the highly non-linear regime where the 
skewness saturates, within our approximation (21), and the 
slope n of the linear power-spectrum shows a very weak de- 
pendence on scale. Therefore, we clearly see in Fig. 9 that 
for smoothing angles 9 S ~ 4' and source redshifts z s ~ 1 
the convergence k s actually probes the intermediate regime 
of gravitational clustering, see also Fig. 6. This sharp fea- 
ture does not show up in the prediction (63) given by the 
spherical-cell approximation because the latter involves an 
integration over redshift along the line of sight. This makes 
the prediction for S^ s smoother and the decline at larger 
angles is shallower since it takes into account the highly non- 
linear scales which are still probed at low redshifts z < z*, 
see also the rise at low z of the typical wavenumber k s shown 
in Fig. 5. This shows that the mean-redshift approximation 
presented in Sect. 4.2 is not sufficiently accurate to obtain 
a good estimate of the skewness S^ s . Of course, at very 
small angles, where the skewness S3 of the density field at 
the typical wavenumber k s (z) only shows a weak variation 
along the line of sight, as discussed above, the mean-redshift 
approximation becomes quite good and it recovers the result 
of eq.(63). 

We may note that at small angles the results of nu- 
merical simulations seem to keep rising rather than to show 
a plateau as for the analytical predictions. This suggests 
that our prescription (21) is only approximate and that 
the skewness keeps slowly increasing within the highly non- 
linear regime. This corresponds to a small deviation from 
the stable-clustering Ansatz. Note that the latter is actually 
a lower bound in the sense that the skewness S3 must either 
remain constant or increase in the highly non-linear regime 
as shown in Valageas (1999) (this only follows from the fact 
that the matter density is positive). An alternative to our 
simple parameterization (21) would be to use a halo model 
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Figure 10. Probability distribution function (pdf) of the 
smoothed convergence k b , V(k s ), for an LCDM cosmology. Var- 
ious panels correspond to various smoothing angles. The sources 
are fixed at a redshift z s = 0.5. Lines correspond to the an- 
alytical predictions as labeled in the Figure: spherical-cell ap- 
proximation (43) (solid line), the mean-redshift approximation 
(47) (dotted line), the log- normal approximation (49) (dot-dash 
line), the stellar-model approximation (61) (dot-long dash line) 
and the Gaussian (dashed line). The predictions obtained from 
the spherical-cell approximation and the stellar-model are almost 
indistinguishable. Data points show the results from numerical 
simulations. 

to evaluate £3 within the non-linear regime. However, we 
shall not investigate such a model in this article since the 
simple parameterization (21) already provides a reasonable 
match to the simulation results over the range of interest, 
as seen in Fig. 9, and, in addition, our prescription seems 
to yield a better match to numerical simulations than the 
halo- model used in Takada & Jain (2002). 

Finally, let us note that higher order moments such as 
the skewness take contributions from the high-K s tail of the 
pdf. Hence they are more affected by the finite size of the 
cutoff than the variance. However we believe that for the 
angular scales considered in our calculations we are not lim- 
ited by the size of the catalogue. We have used several re- 
alizations to probe the underlying mass distribution which 
provides a good handle on the errors introduced sample vari- 
ance. 

6.4 The pdf V(k s ) of the smoothed convergence 

Finally, we investigate here the full probability distribu- 
tion function of the smoothed convergence k 3 . Our ana- 
lytical predictions mainly depend on our prescription for 
the non-linear power-spectrum, taken here from Peacock & 
Dodds (1996), which was specifically tested in Sect. 6.1, on 
our model (21)-(24) for the skewness, which was specifically 
tested in Sect. 6.3, and on our simple parameterization (17)- 
(19) for higher-order cumulants. Therefore, the shape of the 
pdf V(k s ) does not involve new parameters in addition to 
those already tested in the previous sections against numeri- 
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Figure 11. As for the previous Figure, but for source redshift 

Z a = 1. 
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Figure 12. As for the previous Figure, but for source redshift 
z s = 1.5. Notice that although the variance of re s increases with 
source redshift the distribution tends to become more Gaussian 
(see text for discussion). 

cal simulations through the variance and the skewness of the 
convergence. Hence the comparison with simulations of the 
shape of the pdf itself mainly tests the simple form (19) we 
used for the function £(t) (the parameter k being already 
tested through the skewness). 

We compare in Figs. 10- 15 our analytical predictions 
against results from numerical simulations for both LCDM 
and OCDM cosmologies. We consider four smoothing angles, 
8 S — l',2',4',8' and three source redshifts z s — 0.5,1,1.5. 
First, in agreement with Sect. 6.1 we can check that the pdf 
gets broader at higher z s as the variance increases. Second, 
the pdf becomes closer to a Gaussian at higher z s , in agree- 



OCDM 




Figure 13. Probability distribution function (pdf) of the 
smoothed convergence re s , V(k s ), for an OCDM cosmology. Var- 
ious panels correspond to various smoothing angles. The sources 
are fixed at a redshift z s = 0.5. Lines correspond to the an- 
alytical predictions as labeled in the Figure: spherical-cell ap- 
proximation (43) (solid line), the mean-redshift approximation 
(47) (dotted line), the log- normal approximation (49) (dot-dash 
line), the stellar-model approximation (61) (dot-long dash line) 
and the Gaussian (dashed line). The predictions obtained from 
the spherical-cell approximation and the stellar-model are almost 
undistinguishable. Data points show the results from numerical 
simulations. 

ment with the discussion in Sect. 6.1 and with the decrease 
at larger z s of the skewness seen in Sect. 6.3. 

The various line styles correspond to the different ana- 
lytical methods. Thus, we plot the spherical-cell approxima- 
tion (43) (solid line), the mean-redshift approximation (47) 
(dotted line), the log-normal approximation (49) (dot-dash 
line), the stellar-model approximation (61) (dot-long dash 
line) and the Gaussian (dashed line). We can see that the 
Gaussian cannot reproduce the pdf since over these angu- 
lar scales and source redshifts weak-lensing mainly probes 
the non-linear regime so that the pdf is already strongly 
asymmetric. In particular, it shows a sharp cutoff at low 
convergences n s <; « m in an d an extended tail at large posi- 
tive convergences, which follows the shape of the pdf V(8r) 
of the density contrast itself. Next, we note that although 
the log-normal pdf is able to exhibit a large asymmetry and 
provides a significant improvement over the Gaussian, it usu- 
ally fails to reproduce with a reasonable accuracy the results 
from numerical simulations. This discrepancy is especially 
clear at low angles, where the variance and the skewness of 
the density field at scales probed by weak-lensing are large. 
The log-normal approximation fares better at large angles. 
This is due to two effects. First, the variance decreases and 
the pdf gets closer to the Gaussian so that any sensible ap- 
proximation (like the log-normal) which goes to the Gaus- 
sian for small variance would improve in this limit. Second, 
larger angles probe larger scales where the slope n of the 
linear power-spectrum increases and the skewness decreases 
and becomes closer to 3. Then, as noticed in Sect. 4.3 the 
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Figure 14. As for the previous Figure, but for source redshift 

Z a = 1. 
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Figure 15. As for the previous Figure, but for source redshift 
z s = 1.5. 

log-normal approximation is actually very good in the quasi- 
linear regime for n ~ — 1 and 53 ~ 3 (it also coincides with 
our parameterization (17)-(19) in this limit). This explains 
why the log-normal pdf works better at large angles S ~ 8' . 

On the other hand, we can see that our analytical pre- 
dictions based on the simple parameterization described in 
Sect. 3.2 show a reasonably good agreement with numerical 
simulations over all angular scales and source redshifts. The 
improvement over the log-normal pdf is not surprising since 
the model (19)- (24) allows us to follow the dependence on 
time and scale of the skewness 5*3 of the density field for 
any slope n of the linear power-spectrum. This was already 
specifically tested in Sect. 6.3. Then, Figs. 10- 15 show that 
the simple form (19) for the function £(t) which implicitly 



determines the higher-order moments of the density field 
provides a reasonable prescription. In fact, the agreement 
with numerical simulations is surprisingly good in view of 
the simplicity of the model (19)-(24). At large angles our 
prescription seems to overestimate the large-K s tail of the 
pdf. This might be cured by using for the function £(t) the 
exact result derived for the quasi-linear regime (Bernardeau, 
1994, Valageas, 2002), or our simple interpolation (24) may 
overestimate the skewness in the transition regime, although 
this does not seem to be the case for the OCDM cosmology, 
see Fig. 9. However, since the high-K s tail of the pdf for large 
angles 9~>&' may not be of great practical interest we shall 
not try in this article to improve over the model (19)-(24) 
which shows the advantage of a great simplicity. 

From Figs. 10- 15, we note that all approximations 
based on the simple parameterization (19)-(24) yield very 
close predictions. They only show some small differences in 
the far tail of the pdf, at large positive k s , which is very 
sensitive to the details of the model. In particular, the agree- 
ment between the spherical-cell prediction (43) (solid line) 
and the stellar-model prediction (61) (dot-long dash line) 
shows that the spherical-cell approximation (39) is very ac- 
curate and it is sufficient to derive the pdf of the smoothed 
convergence V(ft s ). Therefore, it is not necessary to know 
the detailed behaviour of the many-body correlations £ p : 
the knowledge of the cumulants (<5^) c over spherical cells 
is largely sufficient to predict with a very high accuracy the 
pdf V(k s ) of the smoothed convergence n s . This is an impor- 
tant point since it shows that the measure of the pdf V(k s ) 
would provide a very robust estimate of the pdf V(8r), with 
no degeneracy with the detailed angular dependence of the 
many-body correlations (e.g., whether they are described by 
a minimal tree-model or a stellar-model). 

Eventually, let us note that the numerical computa- 
tion of the pdf from simulation maps suffers from various 
systematics. For a finite size catalogue, clearly the high-K s 
tail cannot continue to infinity. In general the high-K s tail 
shows large fluctuations due to the presence (or the absence) 
of rare overdense (underdense) objects before showing an 
abrupt cutoff. Such effects have been studied in great detail 
for galaxy catalogues (see e.g. Bernardeau et al. 2002) as 
well as for weak lensing surveys (Munshi & Coles 2003). A 
comparison of our numerical results with larger simulation 
volumes will help us to quantify such systematic deviations. 



7 DISCUSSION 

Ongoing weak lensing surveys with wide field CCD imag- 
ing are being used to produce shear maps for areas of order 
10 square degrees. Soon, larger patches, 10° x 10°, will also 
be feasible, e.g., from various surveys including the MEGA- 
CAM camera on the Canada-France-Hawaii Telescope and 
the VLT-Survey Telescope. Such surveys will provide a very 
interesting insight into the dynamics of the universe and the 
clustering in the mass distribution on small angular scales, 
where the density distribution is highly non-linear. They will 
provide opportunities to test our knowledge of gravitational 
clustering at small length scales where no rigorous analytical 
results are presently available. Indeed, unlike galaxy surveys, 
weak lensing surveys provide an unbiased picture of matter 
clustering, as they directly probe the gravitational potential. 
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Results from such future surveys will be especially useful 
where they provide additional data in the form of photo- 
metric redshift information based on the multi-colour data. 
In these cases the additional radial information may assist in 
building up a three-dimensional picture of the dark matter 
structure. 

Statistics of the convergence, k s , can not only be stud- 
ied from the by-product of shear maps generated from weak 
lensing surveys, but can also be studied from the magnifi- 
cation effects of clustered matter which produce variations 
in the image sizes and number density of galaxies across the 
sky (see, e.g., Jain, 2002). Although there are considerable 
difficulties in practice in using these effects, as described at 
length by Bartelmann & Schneider (2001), use of such results 
from wide field surveys can map the large-scale structure and 
help us to quantify its statistics. Recent order-of-magnitude 
analysis of the signal-to-noise ratio as a function of angular 
scale and source redshift have suggested that well-designed 
forthcoming surveys will have high signal-to-noise ratios on 
scales of about 0.1 arcminute to several degrees. This will 
help us to probe the clustering of matter on spatial scales of 
about 50 kpc to 100 Mpc. 

Recent studies of weak lensing have focussed mainly 
on recovering the mass power-spectrum, either from weak 
lensing data alone or from joint analyses of CMB and weak 
lensing surveys (see e.g., Contaldi et al., 2003, for recent es- 
timates). Although such joint analyses can pinpoint the cos- 
mological parameters very effectively, the non-Gaussianities 
we study here form a complementary approach and are able 
to break the degeneracies in estimating the cosmological pa- 
rameters from weak lensing surveys alone. 

Most previous analytical studies of weak lensing statis- 
tics can be divided in two categories. A majority of previ- 
ous works has used a perturbative analysis (Bernardeau et 
al. 1997) which is only applicable in the quasi-linear regime 
and hence will require a large smoothing angle. Although in- 
teresting however, such studies will have limited use for cur- 
rently ongoing surveys as survey areas of order 10 square de- 
grees are necessary to validate such smoothing angles while 
still keeping effects of finite size of the survey area low on 
various statistical quantities. However, given that existing 
CCD cameras typically have diameters of 0.25° — 0.5°, the 
current weak lensing surveys are providing us statistical in- 
formation on small smoothing angles, of order 10' and less. 

Therefore, some earlier studies have also investigated 
the non- linear regime using a hierarchical Ansatz to com- 
pute weak lensing statistics on smaller angular scales (e.g. 
Valageas 2000a,b, Munshi & Jain 2000, Munshi 2000). It 
was shown that the hierarchical ansatz when combined 
with extensions of perturbative calculations such as Hyper- 
extended perturbation theory can provide an accurate pic- 
ture of weak lensing statistics on highly nonlinear scales. 
In this paper we extend these predictions to check their 
validity by comparing analytical predictions with results 
from numerical simulations over a large range of angular 
scales and redshifts. Thus, we have presented two model- 
independent approximations which allow one to obtain the 
pdf of the convergence (and all its moments). The first 
method, the "spherical-cell approximation", simply assumes 
that the relative amplitude of the many-body density corre- 
lations smoothed at a given scale shows a weak dependence 
on the shape of the real-space filter, while the second one, 



the "mean-redshift approximation", makes the further ap- 
proximation that the properties of the density field along 
the line-of-sight can be estimated by considering a typical 
intermediate redshift. Next, we have evaluated the accuracy 
of these approximations by considering a specific "stellar 
model" for the many-body density correlations where exact 
calculations can be performed. We find that in this case the 
spherical-cell approximation is actually very good and can 
barely be distinguished from the exact result. This shows 
that the smoothed convergence k s is a very robust tool to 
measure the statistics of the density fluctuations at a given 
scale and it does not depend on the detailed angular be- 
haviour of the many-body density correlations (it only in- 
volves their overall amplitude at this scale). On the other 
hand, we find that the "mean-redshift approximation" pro- 
vides a reasonable prediction, in agreement with previous 
works (Valageas 2000b), but it does not give a very accu- 
rate prediction for the low-order cumulants like the skew- 
ness at intermediate angular scales, as the properties of the 
density field show significant variations along the line-of- 
sight as it evolves from the linear to the highly non-linear 
regime. A simple model which has been proposed to de- 
scribe the density field from linear to non-linear scales is 
the log-normal approximation. We have applied this model 
to the derivation of the pdf of the convergence based on 
the mean-redshift approximation. Then, we find that while 
the log-normal approximation may be a good description for 
smoothing angles 6 S > 4' it is not sufficiently accurate for 
smaller scales. This is related to the facts that such scales 
are deeper in the non-linear regime and the local slope n 
of the linear power-spectrum P(k) is significantly different 
from — 1, so that the log- normal approximation underesti- 
mates the skewness. Finally, we have applied the three gen- 
eral methods described above (i.e., the spherical-cell, mean- 
redshift and stellar-model approximations) to a very simple 
model for the evolution of the density fluctuations from the 
linear to the highly non-linear regime. We find a good agree- 
ment with the results from numerical simulations for a wide 
range of source redshifts as well as smoothing angles, for 
both LCDM and OCDM cosmologies. 

An alternative approach to model non-linearities in 
gravitational clustering which has emerged in recent years 
is the "Halo model" (see e.g. Cooray A. & Sheth R for a re- 
cent review) . In this approach the dark matter is modeled as 
belonging to halos with a mass function given by the Prcss- 
Schechter formalism and a spatial distribution modeled as 
in Mo & White (1996) or its variants. When supplemented 
by a radial profile for the dark halos such a prescription 
can be used to compute many statistical quantities associ- 
ated with the density field. These models have also been 
applied to model statistics of the convergence field. Our ap- 
proach is complimentary to such methods. It is better-suited 
to the study of weak-gravitational lensing effects which di- 
rectly probe the density field (and its many-body correla- 
tions) and do not involve the decomposition of the universe 
over various classes of objects and virialized halos. In this 
sense, our approach is more natural and the measures ob- 
tained from weak- lensing surveys (e.g., the low-order mo- 
ments of the convergence or the shear) are directly related 
to the density correlations which are the basic ingredients of 
our methods. On the other hand, halo models can provide 
an interesting connection between weak-lensing effects and 
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the properties of virialized halos (so that one might hope 
to derive for instance some constraints on the profiles of 
these halos). However, in our opinion, it is not clear, until 
fully tested, that they can provide a full description of weak- 
lensing effects which also probe low-density regions located 
outside of virialized objects. 

In our studies, both analytical and numerical, we find 
that lower order moments are rather sensitive to the source 
redshift. In this article we only considered the case where all 
sources are located at the same redshift z s but it is straight- 
forward to apply our analytical methods to more realistic 
broad distributions of source redshifts. We plan to inves- 
tigate this point in future studies when relevant numerical 
simulations are available. In agreement with previous results 
we found that while the variance of convergence field (in- 
dependently of smoothing angle) increases with the depth 
of the survey, the distribution of n s tends to become more 
Gaussian, thereby reducing the values of the Sp s parame- 
ters. Nevertheless, the ordering of the Sp B parameters in two 
different cosmologies does not change with source redshifts. 
For any given source redshift the S3 3 parameter is higher in 
the OCDM scenario as compared to the LCDM cosmology. 
The Sp s parameters show a rapid increase at small redshifts 
whereas at large redshifts they vary at a much slower rate. 
Of course, perturbative analysis of the dependence of the 
Sp s parameters on source redshift gives a qualitatively sim- 
ilar dependence. 

Although we have not taken into account the effect of 
the finite size of the catalogue, it is unlikely to have any effect 
on the length scales we have probed in our studies. However, 
the finite size is likely to play an important role in the de- 
termination of the cosmological parameters. It is therefore 
of interest to incorporate such effects in future simulations. 
The effects on lower order moments have been studied an- 
alytically by Munshi & Coles (2003) and we note that the 
higher-order, S p , parameters enter in the expressions for the 
computation of errors for the lower-order moments associ- 
ated with the convergence. 

Noise due to the intrinsic ellipticities of lensed galaxies 
will also need to be modeled. Our analytical results provide 
a strong foundation for such studies and we plan to include 
such features in our future work. The intrinsic ellipticity dis- 
tribution of galaxies can be modeled by assuming a Gaus- 
sian distribution, which depends only on a knowledge of the 
variance in the distribution. To compute the noisy pdf it is 
therefore necessary to convolve the theoretical convergence 
pdf with the noise pdf to produce a resultant containing in- 
formation on both the underlying mass distribution and the 
noise. Initial studies in this direction have been attempted 
by Van Waerbeke (2000) and Jain & Van Waerbeke (1999). 

All analytical calculations regarding weak lensing statis- 
tics rely on the so called Born approximation (Schneider et 
al. 1998). However, its validity at non-linear scales is dif- 
ficult to establish by analytical calculations alone. While a 
weakly clustered dark matter distribution is expected to pro- 
duce small deflections of photon trajectories, it is not clear 
whether the small-angle approximation can be safely used 
at small angular scales which probe large density fluctua- 
tions. However, the agreement of our analytical predictions 
with numerical simulations strongly suggests that the Born 
approximation can be used down to 9 S — 0.4' at least and 
that cases of large deflections of photon trajectories defi- 



nitely have a small statistical significance even on these small 
angular scales. In fact, the theoretical uncertainty is domi- 
nated by the inaccuracy of the simple models used for the 
non-linear evolution of the density field. 

Finally, let us note that our simulations are different 
from more popular ones based on ray-tracing techniques and 
are based on computations of the 3-dimensional shear along 
the lines of sight. The comparison of analytical and numeri- 
cal results for convergence maps are complementary to sim- 
ilar analysis of shear maps presented in Valageas, Barber & 
Munshi (2003). The excellent match between analytical and 
simulation results in both cases increases our faith in both 
the analytical and numerical methods. 
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